Development of anomalous diffusion among crowding proteins 
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In cell membranes, proteins and lipids diffuse in a highly crowded and heterogeneous landscape, where aggregates and dense domains 
of proteins or lipids obstruct the path of diffusing molecules. In general, hindered motion gives rise to anomalous transport, though the 
nature of the onset of this behavior is still under debate and difficult to investigate experimentally. Here, we present a systematic study 
where proteins bound to supported lipid membranes diffuse freely in two dimensions, but are increasingly hindered by the presence of 
other like proteins. In our model system, the surface coverage of the protein avidin on the lipid bilayer is well controlled by varying the 
concentration of biotinylated lipid anchors. Using fluorescence correlation spectroscopy (FCS), we measure the time correlation function 
over long times and convert it to the mean-square displacement of the diffusing proteins. Our approach allows for high precision data and a 
clear distinction between anomalous and normal diffusion. It enables us to investigate the onset of anomalous diffusion, which takes place 
when the area coverage of membrane proteins increases beyond approximately 5%. This transition region exhibits pronounced spatial 
heterogeneities. Increasing the packing fraction further, transport becomes more and more anomalous, manifested in a decrease of the 
exponent of subdiffusion. 
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Introduction 

Anomalous transport of proteins in a crowded environment is best 
discussed in real space, in terms of the mean-square displacement 
(MSD). There, deviations from Fickian diffusion can clearly be 
inferred in a double-logarithmic representation. While computer 
simulations yield the MSD directly, single-particle tracking (SPT) 
is essentially the only experimental method that gives access to this 
quantity G-H Following the trajectories of individual, fiuorescently 
labeled proteins prompts our spatial imagination and facilitates the 
discussion of the underlying mechanisms. Such methods, however, 
require sophisticated and costly image processing, their temporal 
resolution is typically limited to fractions of a second, and long- 
time measurements suffer from bleaching dye due to the intense 
laser light. 

A complementary technique is fluorescence correlation spec- 
troscopy (FCS), which has been established to precisely quantify 
normal protein diffusion at the micro- and millisecond scales.^® 
Compared to SPT, this technique is much simpler and less time 
consuming in terms of performance and data processing, its tem- 
poral resolution is essentially limited by the responsiveness of the 
dye, and the laser power is so low that bleaching is greatly re- 
duced. In particular, FCS offers time windows of several decades 
and is thus a key technique to study anomalous diffusion, e.g., in 
the cytoplasm,!21i2l in cell membranes,^^ and in protein solutions 
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imitating crowding conditions EMS The central quantity in FCS 
is the time-correlation function of the fluctuating intensity of the 
fluorescent light from labeled molecules. A discussion of such a 
correlation function is more subtle than for plain trajectories, and 
it requires certain assumptions on the measurement process itself 
and in many cases even on the nature of the dynamic processes in- 
side the probe. While there is little ambiguity in the case of dilute 
systems showing normal diffusion, the investigation of anomalous 
diffusion is complicated by such prerequisites. Here, we convert 
the FCS autocorrelation function directly to an MSD, from which 
anomalous diffusion can be easily inferred without resorting to a 
particular diffusion model a priori. Further, such an analysis could 
even capture a possible crossover in the MSD from subdiffusion at 
short times to normal diffusion at long times. 

Membrane proteins are naturally exposed to a highly crowded 
landscape, with stark consequences on their transport dynam- 
ics. The cell membrane consists of a heterogeneous lipid bi- 
layer densely packed with integral and peripheral proteins; in ad- 
dition, some of the lipids are organized into microdomains, and 
some of the proteins are tethered to the cytoskeleton. 16 Mea- 
surement of protein transport in cell membranes is thus experi- 
mentally challenging. A pronounced reduction of lateral protein 
diffusivity has been established by numerous experiments using 
SPT,EEZ E] FC S,' 11 ' 12 ' 24 ! 25 ] or fluorescence recovery after photo- 
bleaching.1221 2 ^! There is general agreement that this effect is caused 
by macromolecular crowding .I 3 l 4 l 27 n 29 Many of the experimen- 
tal findings have been interpreted in terms of anomalous diffu- 
sionj 3 ! 1 1 l 12 l 21 ti Hlbut several experiments also display simple Brow- 
nian motion.E 7 l 18 l 22 l It was emphasized that anomalous transport 
should not be considered a universal feature of crowded systems,^ 
and it has been noted elsewhere^ that the transport characteristics 
may depend on the time scale of the observation — which may re- 
solve the controversy. Nevertheless, the relevant crossover time 
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Fig. 1: Experimental system for studying protein crowding on a lipid bi- 
layer membrane. A lipid bilayer of SOPC and biotin-X-DPPE, acting as 
lipid anchors, is formed on a glass substrate; avidin binds irreversibly to 
these anchors. Thus, the surface coverage of avidin is determined by the 
concentration of lipid anchors in the bilayer. A fluorescently-labeled ana- 
log of avidin is included in the avidin mixture for fluorescence correlations 
spectroscopy (FCS). 



scales and, more generally, the source of the anomalous behavior 
have not been identified so far. 

To shed light on the origin of anomalous transport in mem- 
branes, experimental model membrane systems with an adjustable 
degree of crowding are desirable, allowing isolated effects to be 
studied systematically. The hope is to observe both the onset of 
anomalous behavior as well as the crossover to normal diffusion at 
long times as functions of the system parameters. Promising can- 
didates for an experimental realization are synthetic lipid bilayers, 
which offer the advantage that the composition can be controlled 
precisely. Such model membranes have provided insight into how 
diffusion in membranes is affected by immobilized lipids®] and 
lipid phase separation. 31 Yet, the experimental scenario of phys- 
ical molecular crowding in membranes remains relatively under- 
investigated J 27 ! 32 ! 

In the present work, we use FCS to investigate the systematic de- 
velopment of anomalous transport of proteins bound to a crowded 
model membrane on a planar substrate. Specifically, we have stud- 
ied the motion of the protein avidin in supported lipid bilayers as 
function of the protein content (Fig.JTJ. By exploring a wide range 
of avidin surface coverages, this system can be fully character- 
ized in the dilute, intermediate, and crowded regimes, and drastic 
changes in the protein diffusivity are observed. The protein motion 
is clearly inferred from the obtained MSDs, which permit a clear 
distinction between normal and anomalous diffusion. 

Experimental model system 

Our experimental model system consists of a single lipid bilayer 
of SOPC containing a fraction of biotin-X-DPPE, which is sup- 
ported by a transparent glass substrate. 33 The protein avidin irre- 



versibly binds the biotinylated lipid anchors and has confined mo- 
bility in two dimensions on the surface of the bilayer. Avidin has 
two biotin-binding sites that are simultaneously available.^! We 
control the surface coverage of avidin by varying the concentration 
of biotinylated lipid anchors in the bilayer: as the concentration of 
lipid anchors is increased, the surface coverage of avidin increases. 

We take advantage of the non-crystalline nature of avidin at neu- 
tral pH to directly measure how excluded volume affects the pro- 
tein transport. Previous characterization of avidin diffusion using 
fluorescent bleaching techniques showed that at high surface cover- 
ages, avidin does not diffuse on time and length scales comparable 
to the lipids .EE3 In the same lipid bilayer system, the diffusion of 
a lipid tracer below a layer of avidin protein is normal; 33 however, 
the diffusion of avidin itself is difficult to quantify with fluorescent 
bleaching. 

Theory 

Anomalous transport 

The simplest dynamic quantity characterizing the motion of mol- 
ecules is the mean-square displacement (MSD) 6r 2 {t) := (\R(t) - 
R(0)| 2 ) as a function of the time lag t. Simple Fickian diffusion is 
quantified by the diffusion constant D, and the MSD grows linearly 
with time, Sr 2 (t) = ADt for diffusion in a flat membrane. Hindered 
transport is manifested either by a reduction of the diffusion con- 
stant or by a modified shape of the function Sr 2 (t) ; the latter case is 
termed anomalous transport. One generic scenario is a power-law 
relation referred to as subdiffusion, Sr 2 (t) oc f and a < 1, which 
yields a straight line in a double-logarithmic representation. The 
characteristic exponent a is then directly read off as the slope of 
this line. 

Power-law behavior in dynamic correlations is mathematically 
equivalent to the presence of processes with power-law distributed 
characteristic time scalesJSD Generically there is always a slowest 
process, and the subdiffusive behavior is cut off for sufficiently 
long observation times and crosses over to simple diffusion.! 3 ! 37 ! 38 ! 
In the case of cell membranes, a power-law distribution of char- 
acteristic time scales requires a hierarchy of slow processes that 
hinder the molecular motion. One possible source is trapping by 
specific binding J^GU but the origin of the hierarchical ordering is 
not obvious here. Another possibility is that time scales are in- 
herited from a spatially heterogeneous structure, which stems ei- 
ther from the presence of very differently sized macromolecular 
complexes, from self-organized criticality, or which is related to a 
critical phenomenon — as observed for randomly distributed obsta- 
cles.EE] 

The FCS autocorrelation function 

The FCS method consists of autocorrelating the intensity fluctua- 
tions of fluorescent light from a steadily illuminated volume. Since 
the fluorescence is proportional the incident laser intensity W(r) 
and the concentration of labeled molecules c(r, t) fluctuating in 
time and space, the collected light in the detector is determined 
by 

I(t) oc / W(r)c{r,t)d 2 r, (1) 
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where we consider the molecules confined to a plane. The quantity 
measured in an FCS experiment is simply the time-autocorrelation 
function of the fluctuating intensity 6I(t) = /(f) — {/) around the 
mean intensity, 

(SI(t)6I(0)) 
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where the normalization is convention. The intensity fluctuations 
can be expressed as a weighted average of concentration fluctua- 
tions of the labeled molecules. It is straightforward to adapt the 
derivation of Ref. 44 to two-dimensional systems. Then the FCS 
signal 

= ,/d 2 rdVvnr)W(r')S(r-rV) 
(c) [d 2 rtV(r)] 2 
is related to the van-Hove correlation function 



(3) 
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S(r-r',f) = --{6c{r,t)Sc{T' ,0)) . 



(4) 



Its spatial Fourier transform is known as the intermediate scattering 
function^ 



S(q,/) = /s(r-r',/)e 1 i'( r - r V 

and the FCS autocorrelation can be expressed as 
1 Jd 2 g|W(q)| 2 5(q,0 
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where W(q) is the spatial Fourier transform of the intensity profile 
W(r) and 

_ (2nf\W(q = 0)\ 2 _ [Jd 2 rW(r)] 2 

W f<Pq\W{q)\ 2 { > f<PrW(rf ' W 

can be interpreted as the effective number of particles illuminated. 
For a Gaussian profile of the laser beam 



W(r) oc exp ( — 2r /w 



(8) 



where w denotes the beam waist, the filter function iH^q)! 2 °= 
exp(— q 2 w 2 /4) is again Gaussian, and the number of illuminated 
particles is found as N = (c)ttw 2 . 

If only a small fraction of molecules is labeled, the intermediate 
scattering function reduces to the incoherent scattering function 



jq-[R(()-R(0)]\ 



(9) 



For normal diffusion F(q, t) = exp(—Dq 2 t) and again for a Gaus- 
sian beam the FCS autocorrelation then attains the simple form 



G(t) 



1 



1 



Nl+t/T D 
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where tq = w 2 /4D denotes the dwell time. Relaxing the assump- 
tion of diffusion, the incoherent scattering function for two-dimen- 
sional motion is given by F(q,t) = exp [— q 2 Sr 2 (t)/4~\ within a 
Gaussian approximation. Hence, the motion is characterized by the 
mean-square displacement Sr 2 {t) := (AR(/) 2 ), and the FCS auto- 
correlation is then given by 
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A variant of this formula for the three-dimensional situation was 
given during this work by Shusterman et alF^, a generalization to 
polymers is found in Ref. [47] 

The fluorescence conversion in general includes intermediate 
dark triplet states, which modifies the FCS autocorrelation at short 
times. To account for the photo physics, a time-dependent fac- 
tor has been introduced, 48 yet if the triplet lifetime is much shorter 
than the time scales of interest, this factor may be safely ignored. In 
addition, we assume that the fraction of fluorophores in the triplet 
state is sufficiently small and thus negligible for the normalization 
ofG(r). 

Experimental 

Materials 

The lipids l-stearoyl-2-oleoyl-.src-glycero-3-phosphocholine 
(SOPC, Avanti Polar Lipids, Alabaster, Alabama, USA), 
N-((6-(biotinoyl) amino)hexanoyl)- 1 ,2-dihexadecanoyl-.s/?- 

glycero-3-phosphoethanolamine triethylammonium salt 
(biotin-X-DPPE, Invitrogen, Karlsruhe, Germany) and 
2-(4,4-difluoro-5-octyl- 4-bora-3a,4a-diaza-.s-indacene-3- 
pentanoyl)- 1 -hexadecanoyl-jn-glycero-3-phosphocholine 
(/3-BODIPY 500/510C 5 -HPC, Invitrogen) were prepared in 
chloroform. The proteins egg white avidin and fluorescently 
labeled Alexa Fluor 488-avidin (Alexa488-avidin) and the fluo- 
rophore rhodamine 6G were also from Invitrogen. HPLC-grade 
chloroform, acetone, isopropanol, and ethanol were from Carl 
Roth (Karlsruhe, Germany). Reagent-grade NH4OH, 37% HC1, 
and H2O2 were purchased from Sigma (St. Louis, Missouri, 
USA). 

Preparation of the bilayer with bound avidin 

Cover glass slides (24 mmx 24 mm) were cleaned in a 3-stage pro- 
cess and rinsed with deionized (DI) water (MilliQ Corp) between 
each step. First, samples were boiled in a solution of 5:1 volume 
ratio of acetone/DI water for lOmin. Second, samples were boiled 
in 1:1:5 DI water/H202/HCl for 15— 20min, rinsed, then boiled in 
1:1:5 DI water/H202/NH30H. Lipid solutions containing 3 mg to- 
tal lipids were mixed in HPLC-grade chloroform (Roth) in clean 
glass vials. Then, the chloroform was evaporated from the vials 
overnight and then redissolved in HPLC-grade isopropanol to a fi- 
nal concentration of 1.5mg/mL. This solution was used for spin- 
coating by dropping 200 uL of solution onto the clean glass sub- 
strates followed by immediate acceleration of the glass to 3000 
RPM applied for 60s in a Delta 10 spin-coater (BLE Lab Equip- 
ment, Radolfzell, Germany). The lipid-coated substrates were 
dried for 1 hr in a vacuum chamber then placed into Teflon-metal 
chambers that screw together with an O-ring to create a water- 
tight seal to the glass. DI water was added to the chambers to 
hydrate the lipid film. After overnight incubation, the substrates 
were rinsed with DI water and stored at 35— 40 °C for 8hrs, al- 
lowing the bilayer to anneal and to remove any multi-bilayers 



on the substrate J 49 l 50 l The complete removal of possible multilay- 
ered structures was further verified by our previous work using x- 



ray reflectivity and lipid microscopy. 3 ^ The samples were again 
rinsed with DI water, then PBS buffer. Lipid bilayers labeled with 
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0.005mol% /3-BODIPY 500/510C 5 -HPC show normal diffusion, 
even under a fully crowded layer of avidin. Avidin protein solu- 
tions containing 50— 100 ug total protein were added to the sub- 
strates and stored overnight at 35 °C. Before measurement with 
FCS at room temperature, the substrates were rinsed rigorously 
with PBS buffer to remove any non-bound protein. The avidin pro- 
tein solution contained the fluorescent analog of Alexa488-avidin 
with 1:20—1:1000 Alexa488/avidin dye ratio. The dye ratio is op- 
timized at each lipid anchor concentration to maximize the signal- 
to-noise ratio and obtain N ~ 1. 

Fluorescence correlation spectroscopy 

For the FCS measurements, the avidin mixture bound to the mem- 
brane contained a small fraction (0.05—0.001) of fluorescently- 
labeled Alexa488-avidin. ConfoCor2 (Zeiss) with a confocal LSM 
510 microscope was used with a 488-nm Ar laser with power 
165 uW, 6.1 mA. Calibration was done with with Rhodamine 6G 
dissolved in PBS at concentration of lOOnM to determine the di- 
mensions of the focal volume and illuminated membrane area. The 
beam waists obtained varied between w = 0.16um and 0.18 urn. 
All measurements were performed at 21— 22 °C. 

To ensure the z-stability of the sample for long measurement 
times, the drift has to be minimized and controlled. A shift of 
0.4 urn in z-direction results in a 50-fold loss in intensity, and such 
samples cannot be analyzed. During our measurements, the sam- 
ple did not drift more than 0.3 um, and artifacts in the data due to 
drift are not expected. We have also measured the signal from 
an empty chamber with a bilayer but no protein to determine the 
background noise level. The measurements were completed with 
counts per molecule of 1.1 — 1.8, and a noise level < 10 kHz. Mea- 
surements were repeated on at least 2 different samples in at least 
4 different spots. 

Laser scanning microscopy (LSM) 

For fluorescence imaging of the samples, a large pinhole (> 
600 um) was selected for planar imaging of the bilayer. The laser 
power for imaging was 6 mW with pixel time < 2 us. The LSM 
measurements allowed verification that the membranes were intact 
and free of defects. The micrographs in Fig.[7]show averages over 
8 images each. 

Statistical analysis 

To determine whether diffusion is normal or anomalous, a least- 
squares fitting routine using the experimental standard deviation's 
was applied. The parameters N, td, and, in the case of anomalous 
diffusion, tr, were varied, and the autocorrelated data G eX p (t) were 
compared to the diffusion model, G mo( ] e i(?), using the Levenberg- 
Marquardt algorithm in MATLAB. The model, normal or anoma- 
lous, and the corresponding parameters with the smallest x 1 were 
chosen. 

Results and Discussion 

Dilute case 

With the goal of studying the dynamics of avidin on the surface 
of the bilayer over a full range of surface concentrations, we first 
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Fig. 2: Autocorrelated FCS data of the labeled avidin bound to lipid bi- 
layers containing 0.01 (red) and 0.05 (blue) mol% lipid anchors. At such 
low surface coverages, the data are fitted with the normal diffusion model. 
In a separate FCS measurement of labeled lipids (black), normal diffusion 
comparable to the one of avidin is found. The inset magnifies the region 
around the diffusion time. 



examined the dilute case, or the absence of crowding. A very low 
surface coverage was obtained by preparing a bilayer containing 
0.01— 0.05 mol% biotinylated lipid anchors. For such low avidin 
coverage, the measured FCS autocorrelation function G(t) is fit 
well by the normal diffusion model for two dimensions, Eq. l |10| >, 
see Fig. [2] Using D = w 2 /4td, the diffusion constant of avidin 
is D = 2.2 ± 0. 1 um 2 /s and 2.0 ±0.3 um 2 /s at 0.01 and 0.05 mol% 
lipid anchors, respectively. Previous FRAP measurements of the 
homologous protein streptavidin bound to a similar bilayer system 
resulted in comparable values for D. 35 Note that a size-dependence 
of the diffusion constant, such as that described by the Saffman- 
Delbriick diffusion model^l and adaptions thereof is not di- 
rectly relevant here despite the large size difference between lipid 
and protein molecules; avidin is not a membrane inclusion, but it is 
rather bound to the lipid anchors in the bilayer. Using x-ray reflec- 
tivity, we have demonstrated in a previous study 3 ^ that a distinct 
water layer separates the avidin protein layer from the lipid bilayer. 

In a separate control experiment, we have measured the mo- 
tion of the lipids within the bilayer. For this FCS measurement, 
0.005 mol% of lipids were replaced by the fluorophore y3-BODIPY 
500/5 10C 5 -HPC as suggested in Ref. ED We found that the lipids 
in the bilayer with and without bound avidin diffuse normally 
with Diipid = 2.7 ± 0.4um 2 /s. This value is similar to previous 
measurements of the same system using the continuous bleaching 
method. ISlxh us, the bound protein has a diffusion constant compa- 
rable to the lipids in the underlying bilayer, and we conclude that 
at low avidin coverage, avidin transport is essentially limited by 
the diffusion of the lipid anchors. In particular, no signs of anoma- 
lous transport were found. A similar result was obtained in studies 
of diffusion in surfactant bilayers.^ 4 Using a transbilayer peptide 
with streptavidin grafted onto one end, diffusion of the peptide with 
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Fig. 3: Analysis of the FCS autocorrelation function G(t) at a concentra- 
tion of 1 mol% biotinylated lipid anchors. The data are fitted much better 
with an anomalous diffusion model where a = 0.71 (panels b,d) rather than 
with the normal diffusion model (panels a,c). A double-logarithmic rep- 
resentation of the MSD, obtained via Eq. f!4| , yields a straight line with 
slope a (panel f), and is clearly not compatible with normal diffusion, i.e., 
slope 1 (panel e). 



streptavidin is negligibly slowed compared to the bare peptide lack- 
ing streptavidin. 

Anomalous diffusion 

At sufficiently high concentrations of avidin bound to the lipid bi- 
layer, we expect deviations from normal diffusion. Increasing the 
concentration of lipid anchors to 1 mol%, we found that the nor- 
mal diffusion model, Eq. l |10[ >, produces a poor fit as demonstrated 
in Fig. [3^. Instead, the data are fitted much better using a model 
function for anomalous diffusion, 



i(0 = 



N\+{tiT d y 



(12) 



which introduces the exponent a\ here, the fit yields a = 0.71; see 
Fig. [3J). Such a fit implies subdiffusive motion, i.e., a MSD grow- 
ing with a non-integer power in time, Sr 2 (t) oc t a , which follows by 
comparison to Eq. {TTJ derived within a Gaussian approximation. 
An important indicator of the fit quality is to examine the residu- 
als G(t) — G mot iei(f), see Fig. |5J;,d. The residuals for the normal 
diffusion model show a pronounced minimum, which indicates a 
poor fit. 1 1 Examination of the residuals also suggests that the data 
are best fit in the time window 0. 1 ms < t < 300 ms, and that this 
window is limited by statistical noise at long times. 

Let us comment that a superposition of two normally diffusing 
components, 56 



where N — N\ + N2, also allows for a description of the data with 
similarly small residuals. Such a fit, however, is rather susceptible 
to the precise fit range; changing the upper end of the fit window 
from 200ms to 3 s, the ratio Tdll^di varies between 14 and 25, 
while marginally affecting the parameters of the anomalous model, 
Eq. (jT2]». For such large ratios of T d \ /tj2, the FCS correlation func- 
tions resulting from both models can hardly be distinguished in the 
conventional semi-logarithmic representation. Both models, how- 
ever, differ significantly in their long-time decay visible on double 
logarithmic scales. An unambiguous discrimination in the present 
case would require a lower noise level in G(t) of 1%. Nevertheless 
we shall not pursue the two-component fit in the following, since 
the parameters extracted are arbitrary to a large extend and it re- 
mains unclear what the nature of the two independently diffusing 
components should correspond to. 

The usual semi-logarithmic plot of G(t) mainly exhibits the dy- 
namics at the time scale tj, where G(t) ~ 1 /2N. A characterization 
of anomalous diffusion, however, requires sensitivity to a broad 
range of time scales, ideally spanning several decades. In partic- 
ular, the behavior at long times is important, although difficult to 
access due to the relatively small signal of the autocorrelation func- 
tion. Similarly, a least-squares fit is optimized in the vicinity of 
where the statistical noise in G(t) is small.^There have been im- 
portant efforts to improve the quantitative evaluation and analysis 
of autocorrelated FCS data. The maximum entropy method, for 
example, is an adaptation of the^ 2 fitting procedure that considers 
a distribution of characteristic diffusion times.E3E] But as noted 
already, purely subdiffusive transport is equivalent to a continuous 
power-law distribution of characteristic time scales. 

Conversion to the mean-square displacement 

The conventional analysis of an FCS measurement is based on 
model functions that depend on a few parameters such as the diffu- 
sion time tq and the number of fluorophores N in the observation 
volume. Then the diffusion coefficient can be directly extracted 
by means of a fitting algorithm. Such a procedure works as long 
as the model function G mo d e i W reflects the relevant physics of the 
sample under investigation. The selection of an appropriate model 
function is essential for data fitting and the subsequent interpreta- 
tion. If the dynamics of the system under investigation is unknown, 
the proper choice is sometimes unclear and a subjective issue. Usu- 
ally one has to decide a priori whether to use a model for normal 
diffusion, a superposition of diffusing species, or a model for sub- 
diffusion. For the latter, the situation is even more complicated, 
since subdiffusion generically crosses over to normal diffusion for 
sufficiently long times.EMEIl 

Rather than extending existing models by such crossover phe- 
nomena, here we suggest converting the autocorrelated FCS data, 
G(f), directly to the MSD, 



6r 2 (t) 



1 



NG(t) 



- 1 



(14) 



J two-comp\ 



Ni/N N 2 /N \ 

+ . . J , (13) 



1 



N\\+t/T dl l+t/T d2j 



simply by inversion of Eq. {TT}. This eliminates the need to select a 
model function, permitting an unbiased interpretation of the exper- 
iment. Then, a double-logarithmic representation of the MSD re- 
veals the dynamic features at all scales, which can immediately be 
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Fig. 4: FCS autocorrelations and mean-square displacements obtained 
via Eq. (14) for dilute (panels a-c) and crowded (panels d— f) avidin surface 
coverages. As the concentration of crowding protein is increased, the slopes 
of the MSDs decrease in the double-logarithmic representation (panel f). 
Thin black lines indicate fits to the FCS data using model functions for 
normal diffusion (a-c) and subdiffusion (d-f), respectively. 



identified with the naked eye. The beam waist w is accessible from 
independent experiments, and the number N of labeled molecules 
can be inferred from the mean intensity directly, or alternatively, 
from G(t) at short (but not too short) times. Note that the precise 
value of N hardly affects the shape of the resulting MSD at inter- 
mediate and long time scales, and the value of w merely provides 
the conversion to absolute length. 

The inversion of the FCS autocorrelation resolves the tail at long 
measurement times, at the price of magnifying the noise for these 
data. Thus, the experimental challenge is to keep the noise level in 
G(f) low. Statistical noise is minimized by collecting more data: 
longer measurement times give better G(t) statistics. In the present 
work, we take advantage of the stability of the system and have 
measured for relatively long times of approximately 200 s, signif- 
icantly improving the statistics at the time scale of 1 s. An addi- 
tional consideration is the fluorescent signal within the focus area. 
In FCS, the best signal-to-noise ratio for G(t) is obtained by op- 
timizing the concentration of fluorescent molecules such that the 
number of fluorophores diffusing within the measurement volume 
is as close to 1 as possible. For the systematic examination of dif- 
ferent surface coverages, the fraction of dyed protein thus had to 
be optimized for each protein concentration. 

Applying the data inversion method to the example in Fig. [3] 
confirms our earlier conclusion that the avidin diffusion is anoma- 
lous (Fig.[3^-f). The quality of the fit for t > 100 ms can be more 
clearly assessed from visual analysis of the MSD plots compared 
to studying the G{t) plots of the same data (Fig.[5^-b). 
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Fig. 5: The exponent of subdiffusion a (symbols, left axis) decreases as 
a function of the concentration of lipid anchors in the bilayer. The onset 
of anomalous diffusion occurs at a lipid anchor concentration of 0.1 mol%, 
corresponding to an approximately 5% surface coverage of protein on the 
bilayer. Dashed lines, right axis: Upper and lower estimates for the area 
fraction of protein covering the bilayer. 



Development of the exponent of subdiffusion 

To explore the consequences of crowding on the avidin motion, the 
concentration of lipid anchors in the bilayer was systematically in- 
creased. Several independently prepared samples were measured 
in at least four different illumination spots, enabling an estimate 
of the experimental error for the exponent of subdiffusion. In 
the dilute regime with biotinylated lipid anchor concentrations of 
0.01 — 0.05 mol%, the avidin proteins undergo normal diffusion as 
discussed already, and which is further demonstrated by inspection 
of the MSD in Fig. [4]:- 

At a biotinylated lipid concentration of 0. 1 mol%, the FCS au- 
tocorrelation function starts to deviate from normal diffusion; a 
satisfactory fit is obtained using the model for anomalous diffu- 
sion, Eq. (12) , in Fig. [4ji,e. As the biotinylated lipid anchor con- 
centration is increased further, avidin transport on the membrane 
becomes significantly slower, as is apparent from the decreasing 
slope of the MSDs in a double-logarithmic plot in Fig. [4^. Thus, 
the exponent of subdiffusion a decreases monotonically as more 
avidin crowds on the surface of the lipid bilayer (Fig. [3J, and we 
can quantify the development of anomalous diffusion. In the most 
crowded regime, where there is an excess (10mol%) of lipid an- 
chors, we find a = 0.68 ± 0.08. Interestingly, the obtained MSDs 
exhibit single power-law behavior in the accessible time window, 
and the expected crossover to normal diffusion has not yet set in. 
In this case, the conventional method of fitting Eq. \ \ 2\ to the FCS 
function and the conversion to the MSD yield equivalent results for 
the exponent a. 

Thus far, we have described the surface coverage of avidin in 
terms of the concentration of lipid anchors that we used in prepar- 
ing the lipid bilayer. To better depict crowding, the approximate 
area fraction of the membrane covered with avidin may be esti- 
mated based on the geometry of the proteins and lipids. We as- 
sume that each lipid in the bilayer occupies an area of 50 A 2 and 
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Fig. 6: Lateral scanning of a sample with 0.5 mol% lipid anchors yields 
MSDs exhibiting subdiffusion with position-dependent exponents. Inset: 
the exponents as a function of the measurement position. 



each avidin molecule occupies a square of dimension 5 nmx5 nm. 
Depending on whether the avidin molecules are bound by either 
one or two biotinylated lipid anchors, upper and lower estimates of 
the area fraction covered with avidin follow as indicated in Fig. [5] 
Note that the other biotin-binding protein streptavidin was shown 
to bind to a membrane with either one or two biotin anchors. -^The 
onset of anomalous diffusion at 0.1 mol% corresponds to a protein 
area coverage of approximately 3—5%, which is a surprisingly low 
value. 

The transition regime 

By systematically varying the concentration of avidin, we have de- 
termined the onset of anomalous diffusion due to molecular crowd- 
ing. A thorough investigation of these intermediate concentrations 
may elucidate the microscopic origin of the subdiffusive behavior. 
First, we observe that the relative statistical error in the anomalous 
diffusion exponent is significantly higher at intermediate lipid an- 
chor concentrations of 0. 1 and 0.5 mol% (Fig. [5}. Second, evidence 
for heterogeneous structures is found by confocal fluorescence mi- 
crographs in Fig. [7] displaying the distribution of labeled avidin 
bound to the membrane. The lipid bilayer containing 0.5 mol% 
of biotinylated lipid has a more irregular intensity profile in com- 
parison to the dilute and crowded micrographs with 0.05 mol% 
and 10 mol%, respectively. The fluorescence micrographs reveal 
1 — 5 um-sized features in the transition regime. 

These spatial heterogeneities prompted further spatial FCS char- 
acterization of the transition regime. Probing the dynamics on var- 
ious well separated spots allows autocorrelated FCS data locally to 
be collected and local mean-square displacements to be extracted, 
as shown in Fig. [6] All curves are well represented by the subdif- 
fusive law, Eq. {12}, where the exponent interestingly depends on 
the position. Thus our experiments reveal large variations in a on 
a 1 micron scale. Repeated measurement of the same spot after 1 h 
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Fig. 7: Left: fluorescence micrographs at lipid anchor concentrations of 
(a) 0.05 mol%, (b) 0.5 mol%, and (c) 10 mol% reveal spatial heterogeneities 
in the transition regime. Scale bars are 10 pm. Right: intensity profiles of 
the corresponding micrographs, curves are offset for clarity. 



showed no changes in a, corroborating the notion of long-living, 
spatial heterogeneities. 

From the FCS measurements we have collected histograms of 
the local exponents a in Fig. [8] The distributions in the values of a 
suggest qualitatively that there are two classes of transport within 
the transition regime: there are molecules that appear to diffuse 
normally and molecules that show subdiffusive behavior. The his- 
togram for the 10mol% lipid anchor measurements, however, ap- 
pears as a single and very broad distribution. Consequently the spa- 
tial heterogeneities are also present in the crowded regime, where 
at each spot the dynamics is anomalous. Mean and variance of the 
histograms are displayed in the inset of Fig. [8] Treating the two 
classes in the transition regime separately results in significantly 
smaller error bars compared to Fig. [5] 

Conclusions 

We have shown that protein diffusion in crowded membranes dis- 
plays anomalous transport, manifested by subdiffusive behavior of 
the mean-square displacement. The effect of crowding can be stud- 
ied systematically with FCS in a model lipid bilayer membrane 
where the concentration of avidin as crowding agent is well con- 
trolled. Depending on the surface coverage, we observe normal or 
anomalous transport separated by a transition regime. 

We demonstrate an alternative way to analyze the FCS autocor- 
relation function by direct conversion to the MSD. This approach 
is independent of the employed diffusion model and is even capa- 
ble of identifing a possible temporal crossover from subdiffusive 
to diffusive motion. The obtained MSDs, however, exhibit single 
power-law behavior at time scales between 1 ms and 100 ms for all 
concentrations, and we conclude that normal diffusion would take 
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Fig. 8: Histograms of the exponent of subdiffusion a from measurements 
at different positions for lipid anchor concentrations of 0.1 mol% (black), 
0.5 mol% (red), and 10mol% (green). Inset: mean exponents as a function 
of lipid anchor concentration with error bars indicating the variance of the 
histograms. In the transition regime, diffusive and subdiffusive spots were 
treated separately. 



place at longer times beyond the scope of the present set of ex- 
periments. A complementary approach would be to vary the beam 
waist of the laser, enabling a larger area of the bilayer to be mea- 
sured and shifting the sensitivity of the experiment to longer times, 
which could elucidate a crossover mechanism from subdiffusive to 
diffusive motion. 

The transition regime appears to be a mixed phase of micron- 
sized domains of normal diffusion coexisting with domains of sub- 
diffusion. Repeated measurements on the same spot yielded re- 
producible results even after hours, showing that the domains are 
stable and long-living spatial heterogeneities exist on top of the 
membrane; in particular we exclude aging effects. The emergence 
of these spatial heterogeneities correlates with the onset of anoma- 
lous transport. The latter is quantified in terms of the exponent 
of subdiffusion a, which shows a continuous decrease as the de- 
gree of crowding is increased beyond the transition regime. The 
mean-square displacements were extracted directly from the FCS 
autocorrelated data, independent of a priori models for the motion. 
Then different transport mechanisms can be distinguished immedi- 
ately by visual inspection of the MSD. 

A well-established mathematical model for subdiffusion is the 
continuous time random walk (CTRW), where space is homoge- 
neous and a broad distribution of waiting times is assumed from 
the very beginning. In this model, one expects weak ergodicity 
breaking, i.e., time-averaged measurements exhibit aging.^21 Our 
experiments are incompatible with this scenario; rather, they sug- 
gest that the anomalous transport originates from the observed spa- 
tial heterogeneities, similar to the phenomenology of transport on 
percolating systems.^! An analogous conclusion has been made 
recently for a three-dimensional crowded dextran solution.^From 
our experiments we exclude fractional Brownian motion (FBM/^ 
too as a valid description, since it predicts spatially homogeneous 
dynamics. 

One may speculate on the nature of the observed heterogeneous 
structures. The surface oligosaccharide groups of avidin as well 



as its positive surface charge has been postulated to cause lateral 
gel-like network formafionJ^Therefore, jt is possible that protein- 
protein interactions affect the spatial organization of avidin at the 
transition and give rise to the observed heterogeneity. The two- 
dimensional space accessible to membrane proteins may efficiently 
be blocked for topological reasons even at relatively low concen- 
trations of proteins. Thus, the onset of anomalous diffusion being 
at the low surface coverage of 5% supports the idea that protein 
motion is obstructed by spacious, ramified structures. 

In the context of cell membranes, another important question 
is how crowding affects the state of the lipids. Here, we have 
taken advantage of fluid membrane that remains fluid even under a 
very crowded protein layer to eliminate this effect, and anomalous 
diffusion is due to protein crowding alone. The case of integral 
membrane proteins, such as rhodopsin, is more complex,^ and 
further studies are desirable, which introduce protein-induced het- 
erogeneities in the lipid bilayer. 
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